Hall effect in a strong magnetic field: Direct comparisons of 
compressible magnetohydrodynamics and the reduced Hall 
magnetohydrodynamic equations. 

L.N. Martin 1 , P. Dmitruk 1 , D. 0. Gomez 1 ' 2 
1 Departamento de Fisica, Facultad de Ciencias Exactas y Naturales, 
Universidad de Buenos Aires and IFIBA, CONICET, 
Ciudad Universitaria, 1428 Buenos Aires, Argentina. 

2 Instituto de Astronomia y Fisica del Espacio, 
Ciudad Universitaria, 1428 Buenos Aires, Argentina. 
(Dated: September 17, 2010) 

Abstract 

In this work we numerically test a model of Hall magnetohydrodynamics in the presence of a 
strong mean magnetic field, the reduced Hall MHD model (RHMHD) derived by Gomez et al. 
[l|, with the addition of weak compressible effects. The main advantage of this model lies in the 
reduction of computational cost. Nevertheless, up until now the degree of agreement with the 
original Hall MHD system and the range of validity in a regime of turbulence was not established. 
In this work direct numerical simulations of three-dimensional Hall MHD turbulence in the presence 
of a strong mean magnetic field are compared with simulations of the weak compressible RHMHD 
model. The results show that the degree of agreement is very high (when the different assumptions 
of RHMHD, like spectral anisotropy, are satisfied). Nevertheless, when the initial conditions are 
isotropic but the mean magnetic field is maintained strong, the results differ at the beginning but 
asymptotically reach a good agreement at relatively short times. We also found evidence that 
the compressibility still plays a role in the dynamics of these systems, and the weak compressible 
RHMHD model is able to capture these effects. In conclusion the weak compressible RHMHD 
model is a valid approximation of the Hall MHD turbulence in the relevant physical context. 

PACS numbers: 
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I. INTRODUCTION 



The magnetohydrodynamic (one-fluid) model is often regarded as a reasonable description 
of the dynamics of a plasma. One-fluid models are useful in the context of large scale 
dynamics, but when a more detailed description is needed (for instance, when the physical 
context favors the development of small scales) it is most appropriate to consider two-fluid 
models. At this level, the Hall effect, which takes into account the separation between 
electrons and ions, becomes relevant. The Hall length scale distinguishes the range of scales 
at which the Hall effect needs to be taken into account. When the Hall scale is large enough 
to be separated from the dissipation length scale, an approximation known as Hall MHD 
can be used to describe the dynamics. 

The Hall MHD approximation consider two-fluid effects through a generalized Ohm's law 
which includes the Hall current. Whenever one deals with phenomena with characteristic 
length scales comparable or smaller than the ion skin depth c/u{ (c: speed of light; ujc ion 
plasma frequency), the Hall effect cannot be neglected. Among its other manifestations, the 
Hall current causes (in an ideal plasma) the magnetic field to become frozen in the electron 
flow instead of being carried along with the bulk velocity field. Another important feature of 
the ideal Hall-MHD description is the self-consistent presence of parallel (i.e., to the magnetic 
field) electric fields. Hall-MHD has recently been invoked in advancing our understanding of 



phenomena ranging from dynamo mechanisms 



disks 



2J to the physics of turbulent regimes 



magnetic reconnectionj^-5|], and accretion 



In many cases of interest, such as in fusion devices or geophysical and astrophysical plas- 
mas, a strong externally supported magnetic field is present. This external field breaks the 
isotropy of the problem and can be responsible of important changes in the dynamics, such 
as in reconnection regimes or in the development of energy cascades in turbulent systems. 
For one-fluid MHD, the existence of a strong magnetic field is often exploited to yield a 
simpler model: the so-called reduced MHD approximation (RMHD; see refs. 12[ and 13] 
). In this approximation, the fast compressional Alfven mode is eliminated, while the shear 
Alfven and the slow magnetosonic modes are retained [ijj]. The RMHD equations have been 



used to investigate a variety of problem such as current sheet formation 



15 



tionary reconnection 
of turbulence 
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3, 



nonsta- 



18], the dynamics of coronal loops 19|, |20|, or the development 



211 ] . The self-consistency of the RMHD approximation has been analyzed in 
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ref. 22[. Moreover, recent numerical simulations have studied the validity of the RMHD 



equations by directly comparing its predictions with the compressible MHD equations in a 



turbulent regime [23 ] . 



In this paper we test a new model derived in a recent work by Gomez et al. This 
model consists of a system of reduced Hall-MHD (RHMHD) equations derived from the 
incompressible Hall MHD following the same asymptotic procedure, which is employed to 
obtain the conventional RMHD from MHD. The resulting set will describe the slow dynamics 
of a plasma (with Hall currents) embedded in a strong external magnetic field and will 
naturally include new features such as the presence of a parallel electric field. Also we have 
made minor modifications to the model to describe the compressibility effects 24| and we 
tested its performance. 

There are a number of applications for which it is reasonable to consider that the plasma 
is embedded in a strong and uniform magnetic field. Therefore, it is important to find new 
tools and methods to characterize systems where both the Hall effect and a strong mean 
magnetic field are relevant. One of such tools are the direct numerical simulations. The full 
MHD models are computationally demanding and this demand is increased when the Hall 
effect is considered (because of numerical stability issues). On the other hand, the presence 
of a strong mean magnetic field may simplify part of the dynamics. Our goal here is to test 
the RHMHD model 1|, which is aimed at reducing the computational cost. 

The organization of the paper is as follows: Section [II] describes the sets of equations 
for both compressible MHD and reduced Hall MHD, and the codes to numerically integrate 
these equations. In Section HTT1 we present the numerical results: first a comparison between 
incompressible reduced and weakly compressible reduced MHD, second we study the solu- 
tions for different types of initial conditions. Finally, in Section IIVI we list our conclusions. 

II. MODELS AND EQUATIONS 
A. General model 

The compressible Hall MHD equations (dimensionless version) are 



dV Tr 1 J x B (V 2 p*- 1 
V x u> + — T V ' 



dt Ml p \ 2 M|( 7 - 1) 
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V 2 V 1 ,VVT 

P 3 p 



f = -v>n (2) 



^ = VxB- - V0 + 77 V 2 A, (3) 



V • A = (4) 

where V is the velocity field, uj is the vorticity, J is the current, i? the magnetic field, 
p is the density of plasma, A and <fi are the magnetic and electric potential. As shown in 
Eq. ([1]) we used the barotropic law for the fluid, p = cte.p 1 (here p is the pressure), in our 
case we used 7 = 5/3. Ms is the Mach number, Ma is the Alfven number, v and ( are 
the viscosities , rj is the resistivity and e the Hall coefficient. All these numbers are control 
parameters in the numerical simulations. 

We refer to this set of equations as the CMHD equations (Hall compressible MHD). 

B. Reduced models 

The RHMHD model derived in the work by Gomez et al. jlj] is a description of the 
two-fluid plasma dynamics in a strong external magnetic field. The model assumes that the 
normalized magnetic field is of the form (the external field is along e z ) 

B = e g + 5B,\5B\faa<l, (5) 

where a represents the typical tilt of magnetic field lines with respect to the e z direction, 
thus one expects 

Vi«l,9 z «a«l. (6) 
To ensure that the magnetic field B remains divergence free, it is assumed that 

B = e z + V x (ae z + ge x ). (7) 
For the same reason, a solenoidal (i.e. incompressible flow) velocity vector V is proposed, 



V = Vx { V e z + fe x ). (8) 

It is considered that the potentials a(r,t), g(r,t), tp(r,t), and f{r,t) are all of order 
a< 1. Introducing the expressions ((7j) and (jHJ) in the compressible set of equations and 
taking the terms up to first and second order in a the RHMHD model is obtained. 

In a similar fashion, a weakly compressible model can be obtained as well. To this 



end, a slight modification must be introduced to the hypothesis recently described [24J. The 
velocity field, in the more general case, can be decomposed as a superposition of a solenoidal 
part (incompressible flow) plus the gradient of a scalar field (irrotational flow), i.e. 

V = V x (<pe z + fe x ) + W, (9) 

Using the same assumptions as in jl|, but in addition assuming that ip is of order a 2 the 
weak compressible RHMHD equations can be obtained, 

~dt = B °dt + [j ' a] ~ [w ' ^ + uV2lo ' (10) 

^ = B, 9 -^^ + [<p, a] + [b, a] + V V\ (11) 



db „ „ d(u — ej) . 

— = £? /3 P V ^ JJ + b, 6] + p [u, a) + 

-e/3 p [j, a] + /3 p 7/V 2 o, (12) 



— = 5 — + «] - [a, 6] + iA7 2 u, (13) 

where 



<jJ = — V^y?, 

3 = "Via, 
6 = -d v g, 

U = -dyf, 



(14) 
(15) 
(16) 
(17) 
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and the notation [A, B] = d y Ad z B — d y Bd z A is used. 

The only difference between the weakly compressible RHMHD and the original RHMHD 
is in equation (TT2~j) . In the original RHMHD equations (3 P = 1, whereas here (3 P = (3j/ (I+/37), 
with the plasma (3 = 4ttp/Bq. Therefore the compressibility effect here depends on the 
external magnetic field which embeds the plasma. Note that the entire flow compressibility 
is introduced through this parameter, which is a constant of the model. Therefore the 
computational cost is exactly the same as in the original RHMHD model. 

We use this set of equations to study the dynamics of global magnitudes and compare 
them against the results from the full Hall MHD equations and its dependence on the initial 
conditions. 



C. Numerical codes 

We use different pseudospectral codes to solve each set of equations. Periodic boundary 
conditions are assumed in all directions of a cube of side 2irL (where L ~ 1 is the initial 
correlation length of the fluctuations, defined as the length unit). In the codes, Fourier 
components of the fluctuations are evolved in time, starting from a specified set of Fourier 
modes (see following section for the different types of initial conditions), with given total 
energy and random phases. The RHMHD code specifies exactly the same initial fluctuation 
fields (velocity and magnetic field) as prepared for the full MHD code. Initially, parallel 
component fluctuations are set to zero. 

The same resolution is used in both codes in all three directions in most of the runs 
presented here, which have a moderate resolution of 128 3 Fourier modes, allowing many 
different runs to be done with different initial conditions and/or mean magnetic field. The 
kinetic and magnetic Reynolds numbers are defined as R = R m = 1/77, based on unit 
initial rms velocity fluctuation, unit length and non-dimensional values for the viscosity and 
diffusivity. Here we considered R = R m = 400 [y = £ = 1/400, 77 = 1/400) in all of the runs. 
We also considered a Mach number Ms = 1/4, Alfven number Ma = 1, and Hall coeficient 
e = 1/16 in all of the CMHD runs. 



A second-order Runge-Kutta time integration is per 



uated using the standard pseudospectral procedure 



ormed, the nonlinear terms are eval- 



251 ]. The runs are freely evolved for 10 



time units (the initial eddy turnover time is defined in terms of the initial rms velocity fluc- 
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tuation and unit length). Two different values of the magnetic field Bq = 1, 8 are considered 
(in units of the initial rms magnetic fluctuation value). 



III. RESULTS AND DISCUSSION 

We performed simulations with two different initial conditions, one of them maintains the 
anisotropy imposed by the external field while the other does not display any preferential 
direction. 

The first set of initial conditions (anisotropic initial condition) complied with the 
anisotropy imposed by the external magnetic field. This initial condition was generated 
using a two-step process. First, a set of Fourier modes for both magnetic and velocity field 
fluctuations is produced, with amplitudes such that the (omnidirectional) energy spectrum 
is a Kolmogorov spectrum proportional to k~ 5 ^ 3 for 1 < k < 16 and gaussian random phases. 
Modes outside this range in k-space are set to 0. Second, an anisotropic filter is applied, so 
that excited modes for which r n ^(k) > r^k) are initially set to (Here r n ;(k) = l/{kvk) is 
the non-linear time associated with the speed Vk at wavenumber k and ^(k) = l/(k\\VA ) 



is the Alfven time, where Va = -Bq/V^Po is the Alfven speed). This is to ensure that only 



modes which satisfy the RHMHD requirement are excited initially 12Nl4l |26| . Taking the 



Kolmogorov inertial range form Vk ~ k l ^ into account, this condition becomes 

h < C— (18) 

where C is a 0(1) constant. Here we view this inequality as one of the conditions for 
validity of the RMHD approximation. For marginal attainment of the time scale inequality, 
we consider runs with C = 1. The condition means that only large parallel wavelength 
modes (low fcii) are allowed in the initial fluctuations, and becomes increasingly restrictive 
as the value of the mean magnetic field (i.e. Va ) is increased. 

To complete the specification of the initial conditions, the fluctuations are normalized so 
that the initial mean square values of the magnetic and velocity field are both equal to 1 
(unit value). Cross helicity and magnetic helicity are initialized at very small values (remind 
that if the normalized cross helicity is unity there is no turbulent dynamics). Only plane- 
polarized fluctuations are considered in this case, that is, fluctuations are perpendicular to 
the mean magnetic field, e.g., v(x, t = 0) = vj_(x, 0). 
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The runs performed throughout this paper do not contain any magnetic or velocity stirring 
terms, so that the CMHD and RHMHD systems evolve freely. 

In the second case, for the generation of initial conditions (isotropic initial condition), 
we consider initial Fourier modes in a shell in k-space 1 < k < 2 at low wavenumbers, with 
constant amplitudes and random phases. Unlike the previous case, no initial anisotropic 
filter is applied here, so the initial conditions are spectrally isotropic. As in the anisotropic 
case, only plane-polarized fluctuations (transverse to the mean magnetic field) are included, 
so these are (low- to high-frequency) Alfven mode fluctuations and not magnetosonic modes. 

This set of initial conditions represent a more general situation that may arise in an ap- 
plication, where no initial spectral anisotropy is imposed (for instance, they could represent 
better the effect of isotropic driving) . 

The organization of this section is as follows. Subsection IIII Al provides a brief comparison 
of the RHMHD model and its compressible version. In subsection IIII Bl we present the 
numerical results for anisotropic initial conditions. Subsection IIII CI shows the numerical 
results for isotropic initial conditions. 



A. Incompressible RHMHD Vs Compressible RHMHD 

With the introduction of the RHMHD model JJ we continue the line of incompressible 
reduced models (now considering the Hall effect), however it is possible to keep, at very low 
order, the terms responsible of the flow compression. To do this, we slightly changed the 
original RHMHD equations. In this section we compare one model to the other, to see how 
important is the role of compressibility in these systems. The analysis performed on the re- 
sults arising from the different simulations (both anisotropic and isotropic initial conditions) 
show a better fit with the full CMHD model when the compressible version of RHMHD is 
used. Figures ([I]) and (j2J) show different global parameters from the RHMHD simulations 
and its compressible version compared with the same parameters from the CMHD simula- 
tion. Here B = 8, and the initial condition has the anisotropy introduced from the external 
magnetic field. It can be seen that there are slight differences for the incompressible RHMHD 
model whereas a near-perfect fit is observed for the weak compressible RHMHD version. The 
spectra also show noticeable improvements when the compressible version was used. When 
the external magnetic field is changed or using isotropic initial conditions, the weakly com- 
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FIG. 1: Kinetic energy (above) and magnetic energy (below) as function of time for the case Bq = 8 
for anisotropic initial conditions. Column (a) shows the RHMHD model and Column (6) the weak 
compressible RHMHD model. The solid lines corresponds to the CMHD model and the dotted 
lines to RHMHD models. 



pressible RHMHD model still performs better. This result indicates on one hand that the 
relevant parameters of the MHD turbulence in strong magnetic fields are affected by slight 
compressibility effects, and on the other hand that the weakly compressible RHMHD seems 
to fit correctly the above mentioned compressibility effects. 

It is necessary to emphasize that between the RHMHD set of equations and its com- 
pressible version, there is no difference in the computational cost. Therefore, hereafter we 
adopt the compressible version, and refer to it simply with the initials RHMHD. Below we 
compare this (weakly compressible) RHMHD model vs the CMHD model. 



B. Anisotropic initial conditions 

We carry out runs for two different external magnetic field intensities, Bq = 1 and Bo = 8. 
As expected, the global magnitudes studied show that the RHMHD equations fit better when 
the external field is sufficiently intense. Four of these magnitudes (enstrophy, current density, 
kinetic energy and magnetic energy) can be seen in figures |3] and HI 

To obtain a more quantitative measure of the distance between the two solutions 
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FIG. 2: Vorticity (above) and current density (below) as function of time for the case Bq = 8 for 
anisotropic initial conditions. Column (a) shows the RHMHD model and Column (b) the weak 
compressible RHMHD model. The solid lines correspond to the CMHD model and the dotted lines 
to RHMHD models. 




FIG. 3: Vorticity (above) and current density (below) as function of time for anisotropic case. 
Column (a): the results for Bq = 1, the solid lines correspond to the CMHD model and the dotted 
lines to RHMHD model. Column (b): same for the results with Bq = 8. 
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FIG. 4: Kinetic energy (above) and magnetic energy (below) as function of time for anisotropic 
case. Column (a): the results for Bq = 1, the solid lines correspond to the CMHD model and the 
dotted lines to RHMHD model. Column (6): same for the results with Bq = 8. 

(RHMHD and CMHD), we calculate the average distance for the magnitudes shown in 
the figures [3] and HI We define this parameter, for a given magnitude a, as 

J{a RHMHD - a CMHD ) 2 dH 



(19) 



J(a CMHD ) 2 dH 

Table [I] shows this parameter for the magnitudes shown in Fig. [3] and HI Here we can 
see that when increasing the external field the average distance decrease three orders of 
magnitude for the magnetic energy, and two for the kinetic energy. The current density 
decreases two orders of magnitude and the enstrophy one order. There is more than 99,9 
% of agreement for kinetic and magnetic energy when Bq = 8, and approximately 98 % for 
enstrophy and 99,95 % for current density. We can conclude that the RHMHD model is an 
optimal tool for characterize the global magnitudes when the initial condition satisfies the 
anisotropy imposed by the external field. 
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TABLE I: Anisotropic initial conditions 



Magnitude B = 1 B = 8 

£ <V 2> 0.0100301 0.000299110 

£ <B 2> 0.0244613 3.74040 x 10~ 5 

£ <w2> 0.189858 0.0112310 

£ <J2> 0.110552 0.00426486 



In standard RMHD, the parallel components to the external magnetic field do not affect 
the dynamics of the perpendicular fields, since they do not appear in the dynamic equations 
for these fields and in fact act as passive scalars. If the parallel components are initially zero, 
they will remain zero. This situation is, in principle, different in RHMHD, because the Hall 
effect couples the dynamics of the perpendicular components with the parallel components. 
This is evident in equations f fl0|) - ( fl~3]) . If the parallel components are initially zero, they 
may not remain zero at subsequent times, because the Hall term act as a source in the 
evolution equations for those components. Therefore, it is relevant to analyze the behavior 
of the parallel components when starting from initial conditions such that only perpendicular 
fields are not zero. Figure [5] shows the parallel velocity and magnetic field rms (root mean 
square) values. 

The ratio of parallel vs perpendicular components decrease with increasing external field, 
as shown in Fig. [6j In both cases the flow of energy between perpendicular components and 
the parallel component is stabilized. From the Bq = 1 up to Bq = 8 we can see that the 
parallel component goes from about 50 % to less than 6 % for the kinetic energy and 30 % 
to 1.5 % for the magnetic energy. 

Energy spectra for both Bq = 1 and Bq = 8 are shown in Fig. [7] at a time t = 3. We can 
see that when the external field is small, there are slight differences between CMHD and 
the RHMHD model, but these differences disappear when the field increases. For the case 
B = 8, the energy spectra of CMHD and RHMHD are identical. This shows that, when 
the external field is large enough and it maintains the anisotropy imposed by the initial 
conditions, the model reproduces the energy spectrum correctly. Therefore we can say that 
this is a robust model for spectral analysis in a anisotropic case. 
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FIG. 5: Parallel component of the kinetic energy (above) and magnetic energy (below) as functions 
of time for anisotropic case. Column (a) the case Bq = 1 and column (b) the case Bq = 8. The 
solid lines correspond to the CMHD model and the dotted lines to RHMHD model. 
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FIG. 6: Ratio between the parallel component and perpendicular component for kinetic (above) 
and magnetic energy (below). Column (a) the case Bq = 1 and column b the case -Bo = 8. The 
solid lines correspond to the CMHD model and the dotted lines to RHMHD model. 
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FIG. 7: Energy spectrum (kinetic plus magnetic) E(k) for the cases Bq = 1 (above) and Bq = 8 
(below), anisotropic initial conditions, and at t = 3. The solid lines correspond to the CMHD 
model and the dotted lines to the RHMHD model. 

C. Isotropic initial conditions 

In the last subsection we analyze the behavior of the solutions of the RHMHD equations 
whenever the initial condition satisfies the anisotropy imposed by the external field. We 
showed that the parallel fluctuations hardly contributed to the evolution of the system and 
therefore that its evolution is restricted only to the perpendicular plane to the magnetic 
external field. When we consider isotropic initial conditions (and therefore fluctuations 
along the parallel direction) new effects will appear, the most obvious of which are the 
Alfven waves. In this section we aim at quantifying these new effects into the RHMHD 
model and see if this model continues to reproduce correctly the dynamics of the system. 
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FIG. 8: Vorticity (above) and current density (below) as function of time for isotropic case. Column 
(a): the results for Bq = 1, the solid lines correspond to the CMHD model and the dotted lines to 
RHMHD model. Column (6): same for the results with Bq = 8. 

Beside, since only \ow-k modes are excited initially, we expect an energy cascade to develop 
and populate the large wavenumber modes after a few turnover times. We also expect the 
nonlinear activity to be considerably stronger than in the anisotropic case. 

We perform simulations with Bq = 1 and Bq = 8, considering isotropic initial conditions. 
Figure [8] shows the vorticity and current density as a function of time. Here we can see, 
again, a substantial improvement of the fit of the RHMHD curve to the CMHD curve when 
the external field becomes more intense. For Bq = 8 both the overall performance as well 
as the small-scale fluctuations are correctly reproduced, the only difference that appears is 
a slight decrease in the amplitude of the RHMHD curve. However we again emphasize that 
the dynamics is well captured for these global quantities. 

Figure [9] shows the kinetic and magnetic energy as a function of time. Here we see that 
the perfect agreement obtained in the anisotropic case (when B = 8) now gets lost. However 
the difference, as with the vorticity and current, is in an amplitude factor (now, a little bit 
more noticeable). Ignoring the differences due to the amplitude factor, we can see that the 
overall behavior is well described by the RHMHD model, in particular we can see that the 
fluctuations due to Alfven waves are correctly described. 
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FIG. 9: Kinetic energy (above) and magnetic energy (below) as function of time for isotropic case. 
Column (a): the results for Bq = 1, the dark lines correspond to the CMHD model and the light 
lines to RHMHD model (we change the convention here so we can see Alfven waves more clearly). 
Column (b): same for the results for Bq = 8. 

Table [III shows the average distance for global parameters for isotropic initial condi- 
tions. Here we can see that, although the average distances has increased in comparison 
to anisotropic initial conditions, these are maintained below 3 %. Therefore the RHMHD 
model is still very suitable for the description of the system, when the external magnetic 
field is large. 



Magnitude 



TABLE II: Isotropic initial conditions 



B = 1 



B = 8 



£<v 2 > 

£<B 2 > 

£<w 2 > 

£<J 2 > 



0.0409073 
0.0335319 
0.249530 
0.142143 



0.0163489 
0.0244613 
0.0193390 
0.00701709 



In figure [TDl we show the energy spectra (for Bq = 8) at a time in which all the scales 
have been developed. Beside, we show the perpendicular energy spectra. Both spectra show 
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FIG. 10: Energy spectrum (kinetic plus magnetic) E(k) (above) and perpendicular energy spectrum 
E(k±) (below) for the cases Bo = 8 for isotropic initial conditions at t = 3. The solid lines 
correspond to the CMHD model and the dotted lines to the RHMHD model. 

a very good agreement of the RHMHD model with the full CMHD model. 

The energy spectra for the RHMHD model shows a very good fit to the CMHD model, 
validating this model as a useful tool for spectral analysis of plasma embedded in strong 
magnetic fields even if the initial conditions do not maintain the anisotropy of the system. 
The reproducibility of the spectrum is a fundamental feature, because one of the open 
questions is whether the Hall effect might affect the energy inertial range. With this new 
model, which offers significant computational advantages, we expect to be able to address 
unsolved questions like this. 
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IV. CONCLUSIONS 



We have numerically studied the validity of the RHMHD approximation to the full equa- 



tions of compressible 3D MHD (CMHD). As originally proposed in 24J, we took into account 
weak compressibility, without affecting the computational cost. 

The weakly compressible RHMHD model quantitatively reproduces all the global quanti- 
ties when a large external magnetic field is present. The quality of the fit with the 3D results 
show a slight decrease when the initial conditions do not comply with the anisotropy im- 
posed by the external magnetic field. However, this decrease manifests itself in an amplitude 
factor, which does not alter significantly the dynamic evolution of the model. 

We have also observed that parallel component fluctuations to the external magnetic field 
remain small and do not grow with time. 

The energy spectra are always well reproduced, regardless of whether the conditions are 
consistent with the symmetry of the system or not. Since important computational savings 
are offered by the RHMHD model, we can conclude that it can be an excellent tool to study 
the spectral properties of the system, in particular the characterization of the inertial range 
and the dissipation range. 

Acknowledgments 

P. Dmitruk and D. O. Gomez are members of the Carrera Investigador Cientifico CON- 
ICET. Research supported by grants X092/08, X429/08 from University of Buenos Aires, 
PICT 3370/05, PICT 00856/07 from ANPCyT and PIP 11220090100825 from CONICET. 



[1] D. O. Gomez, S. M. Mahajan, and P. Dmitruk, Phys. Plasma 15, 102303 (2008). 

[2] P. D. Mininni, D. O. Gomez, and S. M. Mahajan, Astrophys. J. 584, 1120 (2003). 

[3] F. Mozer, S. Bale, and T. D. Phan, Phys. Rev. Lett. 89, 015002 (2002). 

[4] D. Smith, S. Ghosh, P. Dmitruk, and W. H. Matthaeus, Geophys. Res. Lett. 31, 02805, DOI: 

10.1029/2003GL018689 (2004). 
[5] L. F. Morales, S. Dasso, and D. O. Gomez, J. Geophys. Res. 110, A04204, DOI: 

10.1029/2004JA010675 (2005). 



18 



[6] M. Wardle, Mon. Not. R. Astron. Soc. 303, 239 (1999). 

[7] S. A. Balbus and C. Terquem, Astrophys. J. 552, 235 (2001). 

[8] W. H. Matthaeus, P. Dmitruk, D. Smith, S. Ghosh, and S. Ounghton, Geophys. Res. Lett. 

30, 2104, DOI: 10.29/2003GL017949 (2003). 

[9] P. D. Mininni, D. O. Gomez, and S. M. Mahajan, Astrophys. J. 619, 1019 (2005). 

[10] S. Galtier, J. Plasma Phys. 72, 721 (2006). 

[11] P. Dmitruk and W. H. Matthaeus, Phys. Plasmas 13, 2307 (2006). 

[12] H. R. Strauss, Phys. Fluids 19, 134 (1976) 

[13] D. C. Montgomery, Phys. Scr., T T2/1, 83 (1982) 

[14] G. P. Zank and W. H. Matthaeus, J, Plasma Phys. 48, 85 (1992). 

[15] A. A. van Ballegooijen, Astrophys. J. 311, 1001 (1986). 

[16] D. W. Longcope, and R. N. Sudan, Astrophys. J. 437, 491 (1994). 

[17] D. L. Hendrix and G. van Hoven, Astrophys. J. 467, 887 (1996). 

[18] L. Milano, P. Dmitruk, C. H. Mandrini, D. O. Gomez, and P. Demoulin, Astrophys. J. 521, 
889 (1999). 

[19] D. O. Gomez, and C. Ferro Fontan, Astrophys. J. 394, 662 (1992). 

[20] P. Dmitruk and D. O. Gomez, Astrophys. J. Lett. 527, L63 (1999). 

[21] P. Dmitruk, D. O. Gomez, and W. H. Matthaeus, Phys. Plasmas 10, 3584 (2003). 

[22] S. Oughton, P. Dmitruk, and W. H. Matthaeus, Phys. Plasmas 11, 2214 (2004). 

[23] P. Dmitruk, W. H. Matthaeus, and S. Oughton, Phys. Plasmas 12, 112304 (2005). 

[24] N. H. Bian and D. Tsiklauri, Phys. Plasma 16, 064503 (2009). 

[25] S. Ghosh, m. Hossain, and W. H. Matthaeus, Comput: Phys. Commun. 74, 18 (1993) 

[26] B. B. Kadomtsev and O. P. Pogutse, Sov. Phys. JETP 38, 283 (1974) 



19 



